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AS = 2 and AI = 3/2 Matrix Elements in Quenched QCD * 
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We present preliminary results for Bk, B^ 2 and -Bg /Z from two high-statistics lattice computations. These 
calculations are performed at ft = 6.0 and 6.2 in the quenched approximation, using mean-field-improved 
Sheikholeslami-Wohlert fermionic actions. 
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1. Introduction 



We report on two, high-statistics, quenched lat- 

3/2 3 /2 

tice calculations of Bk, B 7 ' and B s ' . Bk 
parametrizes the K° — A'°-oscillation contribution 
to CP-violation in K decays (e) , leading to a hy- 
perbolic constraint on the summit of the unitarity 

3 /2 

triangle. B 7 ' s are needed to compute the electro- 
penguin contribution to direct CP-violation (e'), 
dominant in the AI = 3/2 channel. Bk and 
Bj g measure deviation from the vacuum satu- 
ration values of the following four-quark, matrix 
elements: 



(K°\0 AS =2\K a ) = -\(0\s lfll 5 d\K°)\ 2 B K , 

with Oas=2 = (srfd)(sj£d), where 7^ L = 1 ± 
75, and in the chiral limit, 

(tt + \0 7 /2 \K + ) -» 7 ](ir + \u 1 5 d\0)(0\s 1 5 u\K+) B 7 /2 

(ir+\0 3 8 /2 \K + ) -> 2(^+|u 7 5 d|0)(0|s7 5 w|iY + )^ /2 , 

where O^ 2 can be written as (sj^d)^^^) + 
(sj^u)(u~f R d) if one forbids penguin contractions 

3/2 

and where 8 is the corresponding color-mixed 
operator. 
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2. Simulation Details 

We describe quarks with the mean-field- 
improved, Sheikholeslami-Wohlert (SW) action 



oSW 



igcsw^ ^2 q(x) PpvVnv q(x) 



with csw = l/ u o an< ^ u o = {|Tr{/ p ;)i and where 
Sp^ is the Wilson fermion action, g, the bare 
gauge coupling, P^v, a lattice definition of the 
field-strength tensor and k, the appropriate quark 
hopping parameter. The parameters of the sim- 
ulation are summarized in Table [j]. We further 
perform full tadpole-improved, KLM rotation of 
quark fields. Though these normalization factors 
cancel in the calculation of B-parameters, they 
specify our renomalization constants. 

3. Operator Matching 

While the matching of quark bilinears is sim- 
ple, the use of Wilson fermions induces mixing 
amongst four-quark operators of different chiral- 
ity To describe this mixing, we use the following 
complete basis of parity-conserving operators: 

r\lat 

r\lat 
^3,4 
Qlat 



'5 



7^ x 7m ± 7^75 x 7m 75 , 
/ x I ± 75 x 75 , 



withTxr= (q'Tq)(q'Tq) lat (a). Then, the match- 
ing can be written, with identifications that will 
be clarified by what follows, 



Oaf= 2 (m) -> Z n 0[ at (a) 
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Table 1 

Simulation parameters. The masses below each k are the corresponding pseudoscalar-meson masses. 

[3 size # cfs. csw K 
6.2 24 3 x 48 188 L442 0.13640 0.13710 0.13745 
a- l (m p ) = 2.h&t\ GeV 780 MeV 570 MeV 438 MeV 
6.0 16 3 x 48 498 1.479 0.13700 0.13810 0.13856 
o~~Vv) = l-96±5 GeV 811 MeV 575 MeV 445 MeV 



and 



oT 
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6 l f{a) 
d[ at (a) 



the chirally subtracted opera- 
; 0[ at (a) + Y.\=iZuO\ a \a) and 



(p) \ / Z 2 2 Z 2 i 

0^(^/2 )^\Z i2 Z Ai 

in terms of 
tors 6[ at {a) e 
d iat {a) _ iat {a) + v _ ^_0^(a), i=2,4. 

We perform matching to the MS-NDR scheme 
at one loop. Since the clover term is O(g), 
we can use the tree-level-clover-action results of 
j]J with modifications appropriate for tadpole- 
improvement and KLM normalization. For the 
coupling, we choose a^ s (fi) defined from the pla- 
quette (|] , identifying the scale /z with that of the 
matching. To estimate the systematic error asso- 
ciated with our procedure, we vary \x in the range 
1/a — > ir/a. 



4. Analysis and Results 

To obtain the desired P-parameters, we con- 
sider ratios of 3-point to two 2-point functions. 
In the limit that the three points are well sepa- 
rated in time, these ratios reduce to: 



Raf=2 



rT 



1 (P(q)\0[ NDR) \P(p)) 



R 
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|(0|pte*|P)|2 



3 (P(q)\Q 2 NDR} \P(p)) 
4 Z^ 5 \(0\P lat \P)\ 2 

l (P(q)\Oi NDR) \P(p)} 
2 Z? \{0\P lat \P)\ 2 ' 



where P is a qq' pseudoscalar meson. We calcu- 
late the ratios for the momenta p — > q = — > 0, 
— » 1, 1 — > lj_ and 1 — f In and for all hopping 
parameter pairs taken from Table U 

To study their chiral behavior, we follow m and 



define the mass and recoil variables 
X 



8 (/j?*) 2 M|, 



3 |(0|P /a '|P)| 2 

Y = P -*X. 

Ml 

We then fit the ratios to 

R(X, Y) = a o + aio X + a 01 Y - 



(1) 



where we neglect both chiral logarithms, which 
are difficult to resolve numerically, and SU (3) f 
breaking terms, which appear to be small for the 
quark masses we consider. 

We find that Raf=2 is well described by a 
linear form in X and Y as shown in Fig. [l] for 
fi = 1/a. At f3 = 6.2 we further find that aoo and 
aio are very small and consistent with zero, as chi- 
ral symmetry requires, and remain so as \i is in- 
creased up to 7T /a. dot should thus give a reliable 
estimate of Bk- At j3 = 6.0, aoo and aio are less 
than 3 and 2a away from zero and smaller than 
the values obtained in Q] with a tree-level SW 
action and boosted, one-loop matching. Taken in 
conjunction with these and other Wilson results 
obtained from less improved actions, our prelim- 
inary findings suggest that discretization errors 
represent an important part of the traditionally 
observed residual chiral violations (please see || 
and d|| for further discussion). 

In the chiral limit, correspond to the 

leading term, aoo, m the expansion of Eq. (Ill) 
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and the study of the F-dependence of R 7 ' 8 is 
less important. Therefore, we consider only the 
p = q = ratios so as not to introduce potential 
momentum-dependent discretization errors. We 
find that the description of the chiral behavior of 
P^g requires a quadratic term in X. 

3/2 

To compare results for Bk and B 7 g obtained 
at different \i and/or /?, we must run them to 
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1 I 1 1 1 1 I 1 1 

= 0.006+0.009-0.008 
= -0.02+0.06-0.08 
= 0.71 + 0.08-0.06 




1 I 1 1 1 1 I 1 1 1 1 I 1 
= -0.019 + 0.008-0.007 
= 0.09 + 0.05-0.06 
= 0.76 + 0.05-0.04 




0.0 0.2 0.4 

Y (p.q) 



0.3 0.4 
Y (p.q) 



Figure 1. Chiral behavior of Raf=2 for fi = 1/a. The different symbols correspond to different p 
different points within a given set correspond to different pseudoscalar meson masses. 



q while the 



Table 2 

B -parameters at 2 GeV in the MS-NDR scheme as a 
function of j3 and the matching scale, /x. 





f 1 


B K 


B 3 7 /2 


„3/2 
n 8 


6.2 


Tr/a 
2/ a 
1/a 


0.7lig 
0.71±| 
0.72+j? 


om±l 

0.58+^ 
0.50t* 


0.8li| 
0.80±| 
0.76±f 


6.0 


Tr/a 
2/ a 
1/a 


0.753 

0.751:4 

0.76± 5 4 


0.55^ 
0.53tj 
0.43 tj 


0.74+;^ 

0.731:3 

0.68^ 



a common reference scale which we take to be 
2 GeV. Running is performed at two loops with 
n/=0. (See Table |.) 

We find that Bk is almost independent of /i in 
the range explored, a-dependence is also found 
to be small though it cannot be excluded given 
the statistical errors and the small deviations of 
a o and a±o from zero at 6.0. Two values of 
the lattice spacing are not sufficient for a proper 
continuum-limit extrapolation and we quote as 
our preliminary result for Bk the [i = 1/a num- 
ber at [3 — 6.2, with the comment that residual 
discretization errors may be small. 

The /^-dependence of B^/ 2 and especially B^ 2 
is significantly stronger than that of Bk- This 
is a result of the rather large matching constants 
and anomalous dimensions. We favor the results 
obtained at the larger values of fx because the 
matching times running is much better behaved 
than at /i = 1/a. Though a-dependence is not, 
on the whole, significant statistically at fixed a/i 



Table 3 

Preliminary results for B -parameters at 2 GeV in the 
MS-NDR scheme (see text). 



B K 


B 3 7 /2 


R 3/2 

B 8 


0.72+1 


0.581S-8 


0.8011H 



or n, discretization errors are difficult to quantify 
because of the large //-dependence. Once again 
we cannot extrapolate to the continuum limit, so 
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the 



we quote as our preliminary results for B 7 g 
/i = 2/ a numbers at (3 — 6.2 to which we add 
systematic errors to account for the observed /it- 
dependence. 

Our results are summarized in Table [|. Fur- 
ther discussion of systematic uncertainties must 
be postponed for lack of space. For further com- 
parisons with other recent results, please see||. 
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